function R1=findR(Q1,I,dt,lambda,n,fL1,FL1)
R1=zeros(I,n+1);%加一列给k=0的情况
%Yi-1的理论最大值是前i-1列的病人数的和,我们把范围放大到n，这并不影响结果。
%m的范围为X的m的范围,我们与Q1保持一致把范围放大到n，这并不影响结果。

%k=0时的概率，即求R1(i,1)=R(i,0),求R1的第一列
%发现利润与R的第一行没关系，也许可以直接不算第一列，同时放弃R给第一列的设定。
%求R1的第1行
%计算R1(1,1)
for m=0:n 
R1(1,1)=R1(1,1)+(1-FL1(m+1))*Q1(1,m+1);
end
%计算R(1,2:n+1)
for k=1:n %k的理论最大值为n，即所有病人都在i时期及之前来了，且未得到治疗。
  for m=0:n        
    if m-k+1>0
      R1(1,k+1)=R1(1,k+1)+fL1(m-k+1)*Q1(1,m+1);
    end
  end
end

%求R1的第2到n行
for i=2:I
  %计算i行第1列
  for k2=0:n%此处k2指k~
    for m=0:n   
      R1(i,1)=R1(i,1)+(1-FL1(m+k2+1))*Q1(i,m+1)*R1(i-1,k2+1);
    end
  end
  
  %计算i行第2到n+1列
  for k=1:n
    for m=0:n   
      for k2=0:n%注意这里的k2指k~
        if m+k2-k+1>0   
        R1(i,k+1)=R1(i,k+1)+fL1(m+k2-k+1)*Q1(i,m+1)*R1(i-1,k2+1);
        end
      end
    end
  end
end


   

    
